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Abstract 



The standard numerical approach to determining matrix elements of local opera- 
tors and width of resonances uses the finite volume dependence of energy levels and 
matrix elements. Finite size corrections that decay exponentially in the volume are 
usually neglected or taken into account using perturbation expansion in effective field 
(D I theory. Using two-dimensional sine-Gordon field theory as "toy model" it is shown 

that some exponential finite size effects could be much larger than previously thought, 
potentially spoiling the determination of matrix elements in frameworks such as lattice 
QCD. The particular class of finite size corrections considered here are //-terms aris- 
ing from bound state poles in the scattering amplitudes. In sine-Gordon model, these 
can be explicitly evaluated and shown to explain the observed discrepancies to high 
precision. It is argued that the effects observed are not special to the two-dimensional 
setting, but rather depend on general field theoretic features that are common with 
models relevant for particle physics. It is important to understand these finite size 
corrections as they present a potentially dangerous source of systematic errors for the 
determination of matrix elements and resonance widths. 



H 1 Introduction 



The matrix elements of local operators (form factors) are central objects in quantum field 
theory. Indeed they are the subject of considerable interest in lattice QCD, where - among 
other applications - they are relevant in describing the weak decays of hadrons. Due to the 
Maiani- Testa no-go theorem [1], it is necessary to extract them using finite-size methods, 
which are the subject of a seminal paper by Lellouch and Liischer [2] (see also [3]). A closely 
related problem is description of resonances, whose decay width can also be extracted from 
finite-size data, for which the first proposal was made by Liischer [4]. 

Two-dimensional integrable quantum field theories provide an ideal testing ground for 
these ideas because in many such models the exact analytic expressions of the form factors 
are known. First, the S matrix can be obtained exactly in the framework of factorized 
scattering developed in [5] (for a later review see [6]). It was shown in [7] that it is possible 
to obtain a set of equations satisfied by the form factors using the exactly known scattering 
amplitudes as input. The complete system of form factor equations, which provides the 
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basis for a programmatic approach (the so-called form factor bootstrap) was proposed in 
[8]. For a detailed and thorough exposition of the subject we refer to [9]. 

On the other hand, for two-dimensional theories there is a very efficient alternative to 
lattice field theory to evaluate finite size spectra and matrix elements called the truncated 
conformal space approach, which was first proposed by Yurov and Zamolodchikov [10]. 

Recently we used this framework to perform a detailed analysis of resonances and local 
operator matrix elements. Resonances were studied in [11], where it was argued that they 
can be extracted with a much better accuracy from level splittings than from the slope 
method originally proposed in [4]. It was realized that the proposed method (dubbed the 
"improved mini-Hamiltonian" method) is essentially equivalent to determining the matrix 
element of the interaction term responsible for the decay, and in [12, 13] we explored an 
approach to treat general matrix elements of local operators, which is an extension of 
the Lellouch-Liischer approach. A proof given in [12] shows that results obtained in this 
framework are valid to all orders in the inverse volume 1/L, i.e. up to corrections decaying 
exponentially with the volume. We also realized that these so-called "residual" finite size 
corrections could play an important role, but the tools to study a particular class of these, 
the so-called /x-terms were only developed later by Pozsgay in [14]. 

In this work the issue of yU-terms is addressed in more details. The testing ground is 
a specific version of the well-known sine-Gordon theory, the so-called /c-folded model [15], 
mainly because all the necessary techniques are well-developed and the model is thoroughly 
understood from the theoretical point of view. This work can also be considered as a 
refinement of the form factor study performed in [16], where the relevance of //-terms was 
pointed out, but their detailed analysis was omitted. 

The main goal of this paper is to present the issues in a way that leads to conclusions 
which are expected to be relevant to a wide class of models, and in particular to the ongoing 
efforts to determine resonance parameters and decay matrix elements in lattice QCD. 

The outline of the paper is as follows. In section 2 the necessary details concerning the 
two-folded sine-Gordon model are introduced. Section 3 presents the formalism of finite 
size form factors, first the description valid to all orders in 1/L and then determining the 
leading exponential corrections a.k.a. the /i-terms. The theoretical results presented here 
are compared to numerical data extracted from TCSA in section 4. Section 5 discusses 
the relevance of these findings to the description of resonances, and in section 6 the con- 
clusions are formulated. The detailed and rather bulky formulae for the exact form factors 
(determined from the bootstrap) are summarized in Appendix A. 



2 The two-folded sine-Gordon model in finite volume 

The classical action of sine-Gordon theory is 

^ = y"rf2x Qa^$a^$ + ^cos/3$^ (2.1) 

A useful representation of the model at quantum level is to consider it as the conformal 
field theory of a free massless boson perturbed by a relevant operator. In this framework, 
the Hamiltonian can be written as 

H = J dx^: {dt^f + {d^^f : -X J dx : cos /3$ : (2.2) 

where the semicolon denotes normal ordering in terms of the modes of the A = massless 
field. Due to anomalous dimension of the normal ordered cosine operator, the coefficient 
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A has dimension 

A ~ fmass 



2-/32/47r 



and it sets the mass scale of the model. The genuine coupling constant is /3 which for later 
convenience is reparametrized as 

^ Stt - /32 

The spectrum consists of a doublet of solitons and their bound states, which are called 
breathers. The scattering amplitudes of this model are briefly reviewed in Appendix A.l. 
As usual in two-dimensional kinematics, the on-shell energy-momentum two- vector {E,p) 
of a particle with mass m is parametrized with the rapidity 6: 

E = m cosh 6 , p = m sinh 6 

Lorentz boosts correspond to shifting all rapidities by the same constant. 

In finite volume, it is possible to choose quasi-periodic boundary conditions for the field 
leading to a set of possible local theories labeled by an integer k called the folding number 
[15]. To define the k -folded theory SG{f3, k) we take the sine-Gordon field $ as an angular 
variable with the period 

271 

$ ~ $ + — -A; (2.3) 

When the space is a finite circle with circumference (in this case also the volume) L, i.e. 
X ^ X + L, this implies the following quasi-periodic boundary condition for the field 

2n 

^{x + L,t) = ^{x,t) + —km , m e Z (2.4) 

The classical ground states are easily obtained: 

$ = ^n, n = 0,...,A;-l (2.5) 

which shows that the condition (2.3) corresponds to identifying the minima of the co- 
sine potential with a period k . In the infinite volume (L = oo ) quantum theory these 
correspond to vacuum states \n) which have the property 

2tt 

{n\<^{x,t)\n) = —n n = 0,1, . . . ,k - 1 (2.6) 

/-^ 

These states are all degenerate in the classical theory and also at quantum level when 
L = oo ; however, tunneling lifts the degeneracy in finite volume L < oo . The spectrum of 
breather multi-particle states is obtained in k copies corresponding to the k vacua, and in 
finite volume tunneling splits the degeneracy between these states by an amount of order 
e~*^^ where M is the soliton mass. In the /c-folded model, the local operators of most 
interest are the exponential fields 

In the sequel the value of k is set to 2 so the model considered is the 2-folded sine-Gordon 
theory SG{l3, 2) . 
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3 Form factors in finite volume 



3.1 Corrections to all order in 1/L 

In this work we restrict our attention to states containing only breathers (i.e. no sohtons). 
As breathers are singlets and therefore scatter diagonally, the formulae derived by Pozsgay 
and Takacs in [12, 13] are directly applicable. The first ingredient is to describe the multi- 
breather energy levels corresponding to the states 

whose finite volume counter part we are going to label 

. . . , /Ar})ri...rjv,L 

where the Ik are the momentum quantum numbers. In a finite volume L, momentum 
quantization is governed (up to corrections decreasing exponentially with L) by the Bethe- 
Yang equations: 

Qk{Oi, ...,9n)= rrir^L sinh 6*^ + ^ 5.r,rfe (^fc - ^j) = 27r4 4 G Z (3.1) 
where the phase-shift is defined as 

S,,{d) = e'^'^'^^^ 

These equations are nothing else than an appropriate extension of the description of two- 
particle scattering states in finite volume [17, 18]. In general quantum field theories this 
description is only valid under the inelastic threshold, but due to integrability it can be 
extended to all multi-particle states regardless of their energy. 

In practical calculations, because breathers of the same species satisfy an effective 
exclusion rule due to 

^^^(0) = -1 

it is best to redefine phase-shifts corresponding to them by extracting a minus sign: 

SO that 6sr{0) = can be taken for all s,r and all the phase-shifts can be defined as 
continuous functions over the whole real 6 axis. This entails shifting appropriate quantum 
numbers Ik to half-integer values. Given a solution ^i, . . . , ^tv to the quantization relations 
(3.1) the energy and the momentum of the state can be written as 

N 



E = m^j. cosh 9k 



k=l 

^ ~ 2ix 

P = ^ m^^ sinh 9k = — 

k=l k 

(using that - with our choice of the phase-shift functions - unitarity entails Ssr{9) + 
^rs{—9) = 0). The rapidity-space density of n-particle states can be calculated as 

Pr,...rM,---,On)=detJ^-^ , jin) ^ 9Qk{9u^. . , 9^) ^ ^ _ _ _ ^ ^ (3,2) 
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In infinite volume the matrix elements of local operators between multi-particle states 
composed of breathers can be expressed in terms of the elementary form factor functions 



...,0^) = {0\O\BrM . . . B^M) (3.3) 

using crossing symmetry 

+ disconnected contributions (3.4) 

where disconnected contributions only arise when there are particles in the two states that 
have identical quantum numbers and momenta. 

In terms of the elementary form factors (3.3), the finite volume matrix elements of local 
operators between multi-particle states can be written as [12] 



S1...SM iWi^ • • • ' 4/}|C(0, 0) . . . , In}) 



\J Pri...rjv(6'l, . . . , ON)Psi...SMi(^'l^ • • • > ^m) 



0(e-^^) (3.5) 



which is valid provided there are no disconnected terms. This is essentially the extension 
of the Lellouch-Liischer formalism [2] to general multi-particle states. The finite volume 
quantization relations in (3.1) and the expression of density factors p in (3.2) are specific 
to the setting of two-dimensional integrable field theories, but the relation (3.5) which 
states that finite volume matrix elements differ from their infinite volume counterparts by 
normalization factors given by the square roots of densities of states is generally valid in 
any quantum field theory (with a nonzero mass gap), as can be seen by considering the 
way this expression was derived in [12]. 

Note the absolute values in (3.5) that are necessary to take into account that the finite 
volume and infinite volume phase conventions differ. In our framework it is known how to 
compensate for the difference. The finite volume eigenvectors obtained from TCSA can be 
chosen real, as the Hamiltonian is a real symmetric matrix, and then all matrix elements 
of exponential operators turn out to be real. On the other hand, apart from the i factors in 
(A. 3) the only complex function is the minimal two-particle form factor f^{0), whose phase 
at a real value of the argument is equal to the square root of the S'-matrix. Therefore the 
above relation can be rewritten as 



n...sM iUi^ • • • ' ^mIIC'IO, 0) . . . , In}) r-i...ri^,L 

X (3.6) 



FZ...s,r,...rA^'M + ^Vr, . . . , + ZTT, ^1, . . . , ^^) 



\/ Pri...rjv(6'l, . . . , ON)Psi...SMi^l^ ' ' ' ' ^A/) 
l<k<l<M l<k<l<N 

where the only remaining ambiguity is a sign corresponding to the choice of the branch of 
the square root (note that this ambiguity is also manifest in the TCSA eigenvectors since 
demanding their reality does not fix their sign). Note also that the values 

are real whenever 6 is real. 
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Disconnected terms may appear when there are two breathers of the same species in 
the two states whose rapidities exactly coincide. This can happen in two cases: when the 
two states are identical, or when there is a particle with exactly zero momentum in both 
of them [13]. Here we only quote the latter since that is needed for our computations. 
Defining the function 

J^k,ii6i, • • • , ^fel^i, ■ ■■ = 

lim F 11 1 (in + 9[ + e, . . . , in + 9i + e, in — 9i + e, . . . , in — 9[ + e, iir + e, 

2fe+2i+2 

0,9u...,9i,-9i,...,-9i) (3.7) 
one can write 

n^({/(, ...,I'„0, -II -/ai<^'|{/i, 0, -Ii, -/i}) n^,L (3.8) 

2fc+l 21+1 

. . . , 9',\9,, mL Fn...i (^vr + 9[,...,tn- 9[, 9^, 



X 



P2/t+l(^l; • • • ; ^fc' 0' • • • ' -^l)P2Z+l(^l: • • • , 6';, 0, -9^ • • • , -9i) 

(phase factor) + 0(e~''^) 

where the compensating phase factor is entirely analogous to the one in (3.6). 

The above predictions for finite volume energy levels and matrix elements are expected 
to be exact to all (finite) orders in 1/L [12, 13] (note that the exponential corrections are 
non-analytic in this variable). 

3.2 /i-terms: the leading exponential corrections 
3.2.1 yU-terms for energy levels 

Using the ideas in [14], we can model a finite volume B2 state as a pair of Bi particles with 
complex conjugate rapidities 

B2{9) ^ Bi{9 + iu)Bi{9 - lu) 

where u can be obtained by solving the -Bi-Bi Bethe-Yang equations (written here in ex- 
ponential form): 

gimiLsinh(e±»)^^^^^2m) = 1 (3.9) 

where 

sinh 6* 2 sin TT^ 
smh 9 — t sm n^ 
The solution for u has the large volume behavior 

M ^ !l + tan ^e-^^^lLcoshe ^3_^q^ 

with 



2 I 2 ^2 ■ < 

f^n = ymf- — = mi sm — 

Therefore, for L — )■ 00 one obtains the usual bootstrap identification 

B2{9) ^Bi{9 + ini/2)Bi{9 - ini/2) 
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(a) B2 



(b) B, 

Figure 3.1: /i-term diagrams 



and the momentum and energy of the state tends to 

mi sinh(6' + m) + mi sinh(6' — m) = 2mi cosnsinh^ — )■ m2sinh6' 
mi cosh(6' + m) + mi cosh(^ — iw) = 2mi cosM cosh^ — m2Cosh6' 

where 

m2 = 2mi cos — 

is the mass of B2 in terms of that of i?i, consistent with (A. 2). For a stationary i?2, the 
energy dependence on the volume is 

E2{L) = 2miCOSM ~ m2 - (7n)Vne~''"^ + • • • (3.11) 

where the second expression shows the asymptotic behavior vahd for large enough L. One 
can also write 



/ill — ''^1 sinwjg 



where the -B1-B2 Bi fusion angle is 



= ( 1 - I 



The asymptotic behaviour in (3.11) coincides with that predicted by Liischer's finite 
volume mass formula [19], which was worked out in more detail for 2-dimensional quantum 
field theory in [20] and is illustrated graphically in Fig. 3.1 (a). Here 7^1 is the three- 
particle coupling defined in (A. 8). The kinematical variables can be represented by the 
mass triangle shown in Fig. 3.2, with /i^^, being the height of the triangle perpendicular to 
the side mc. Note that can be much smaller than other mass scales in the theory if the 
binding energy is low, i.e. when ttlc ~ m^, + m?,. 

A stationary third breather B-^ can be modeled as a composite of three Bi particles 

B^i{e = 0) ~ 5i(m)fii(0)fii(-m) 

where the (nontrivial) Bethe-Yang equations read 

gimiLsinh(±m)^^^^^^^)^^^(^2m) = 1 (3.12) 

leading to the asymptotic expression 

^ ^ 4sin< + 2tanvre 

^ 2cos<-l ^ ' 
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Figure 3.2: Mass triangle, fusion angles and /i parameter 



The energy of the level then becomes 

EsiL) = mi(l + 2cosm) ~ - 2 (tL) V^e"'^'^'' 
where the three-particle coupling 7^2 is defined in (A. 8) and 



(3.14) 



1A2 — sin = 1712 sin 



with the fusion angles 



u 
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^(1-0 



U23 = TT { 1 - I 



,1 



The asymptotic behavior in (3.14) is again consistent with the expression given by Liischer 
[19, 20], the factor of 2 is due to the presence of two processes contributing the same 
amount as shown in Fig. 3.1 (b). 

3.2.2 yU-term corrections for the form factors 

Using the results of [14], these can be computed by continuing (3.6) to complex rapidities 
that correspond to representing B2 and B3 as composites made out of Bi. One obtains for 
vacuum- i?2 matrix elements the form 



(0|(9(0,0)|{/})2,, 



lU] 



(3.15) 



Pii{0 + iu,6 — iu) 
where 6', u is the solution of (3.9) with the correct momentum, i.e. 

miL sinh(6' + iu) + rriiL sinh(6' — iu) = 21x1 

and the dots denote further (and generally much smaller) exponential corrections; the ± 
sign accounts for the remaining phase ambiguity from the square root. One can similarly 
evaluate Bi — B2 and BiBi — B2 matrix elements using 
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_{{I'}\O{0,0)\{I})2,i 



{{I[,I',}\OiO,0)\{I})2,i 



. ^/Su{2^u)Ff^^{i^T + 6', 6- m, 6 + iu) 
rriiL cosh 9' pii{d + iu,6 ~ iu) 



(3.16) 



± 



S^^{2^u)S^^{9[-9',] 



Pii{0[,9'2) pu{9 + iu,9 -iu) 



X 



^iiii(^^ + 92,i7r + 9[,9 - iu,9 + iu) 
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For large enough L, we can use (3.10) and keep the leading term of these expressions in 
the limit -u — > 7r^/2 to obtain 



(0|0(0,0)|{/})2,L 
,{{I'}\O{0,0)\{I})2,L 

u{{i[,i^}\o{o,o)mhL 



7712 L cosh ^ ^ 



(3.17) 



±i 



\/ TJiiL cosh 9' cosh 6 



+ O (^e~^?i^ 



Si,{e[-e',)Fg{z7r + e',,z7r + e[,i 



X 



Puid'i - 92)m2Lcoshe 
+0 (|e-^"^) 

where we also used (A. 7). This is just what we obtain from the general expression (3.6). 
In performing this calculation, the singular behavior of the factor S'ii(2m) for u — )■ 71^/2 
is canceled by a singular term from the denominator factor 

,x _x m2-Lcosh6' 
Piiiu + iu,0 — lu) = — ; T— + regular terms 

For matrix elements involving a stationary we obtain 



(0|0(0,0)|{0})3,L 

,{{i'}\o{0M{0}h,L 

n({/(,/aiC?(0,0)|{0})3,L 



_^ . ^ySn{2iu)Sll{iu)Ff^^{-iu, 0, +iu) 

A/piii(+m,0, -iu) 
_^ ^ySn{2iu)Sn {iu)F^^^ {in + 9', -iu, 0, +w) 

rriiL cosh 9' pui{+iu, 0, —iu) 



(3.18) 



+ 



±1 



^iiiii(^7r + 6*2, in + 9[, -iu, 0, +iu) 



Sni2tu)Siii9[-9'2)Si,itu) 
pn{9[, 9'^) piii(+m, 0, -iu) 



+ ... 



where u solves (3.12) with the asymptotics (3.13). There is an exception to the second 
formula, however: when /' = 0, the Bi is stationary and there is a disconnected piece, 
which according to (3.8) leads to the modification 



_({0}|O(0,0)|{0})3,, 



where 



a/ 5*11 (2m) 5*11 (m) J'fJ'3(-m) + miL Sii{2iu) F[1{ iu, -\-iu) 



rriiL cosh 9' pm {+iu, 0, —iu) 



+ . . . 



(3.19) 



J^^,,{9)=\unF^,,,{zn + e,9,0, 



Using the exchange property 



^V...Vfc+i..i„ (^1' • • • ' ^k, 9k+l, ...,9n) 



Si 



ikik+1 



6'fc+l)-^f..ife+iifc..i„(6'l, • • • ,9k+l,9k, . . .,9n) 



(3.20) 
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we can write 

'-0 (a\ CI h\ -nO 



KM = si-e)K3id) 

where ^{6) = hm F^^^-^{i7r + e, 0, 6, 



which explains the presence of the factor 5'ii(m) in the coefficient of J-'fg in (3.19), in 
contrast to formula (3.8). (The factor 1^5*11 (2m), common to both terms results from 
the analytic continuation of the compensating phase factors discussed in the previous 
subsection, as all similar factors in other cases do as well.) 

It can again be verified that for L large enough, taking the limit u ^ again results 
in agreement with what is expected from (3.6): 

(0|0(0,0)|{0})3,L = ±:^i^ + O (e-^-^) (3.21) 



({/'}|O(0,0)|{0})3,^ = ±z^^Ei£l£L + O (e-^?^^) 



n(W,/^}|O(0,0)|{0})3,^ 



a/ niiL cosh 9' m^L 



+0 (e 



X 



Note that the additional term in (3.19) drops out because the factor S'ii(2m) has no pole 
to compensate for the singular behavior of the denominator which is of the form 

Piii(+m, 0, —iu) ~ — — r + less singular terms 

[u — nt,)^ 



3.3 Beyond the leading ^u-terms 

It is important to consider subleading exponential corrections. First of all, the breathers 
also have /x-terms corresponding to their description as soliton-antisoliton bound states, 
with the characteristic scale 

In addition to the /i-terms there are corrections which correspond to particle loops winding 
around the space-time cylinder. They are called F-terms because they depend on the 
forward scattering amplitude F continued to rapidity difference with imaginary part in/ 2 
[19] and are suppressed by the mass of the particle propagating in the loop. Last, but not 
least, there are also the tunneling corrections related to the presence of multiple vacua {k = 
2 in the case considered below), suppressed by the soliton mass M, which are essentially 
vacuum versions of F-terms. These contributions were evaluated using the dilute instanton 
gas approximation in [15]. F-terms contributions are illustrated in Fig. and their scales 
are compared explicitly to //-term scales at the beginning of the next section. 



4 Determining matrix elements of local operators: the 
effect of yU-terms 

To evaluate the form factors numerically, we use the truncated conformal space approach 
(TCSA) pioneered by Yurov and Zamolodchikov [10]. The essential idea is to truncate 
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(a) (b) 

Figure 3.3: (a) F-term correction to the mass, the filled circle denotes the forward scattering 
amplitude evaluated at ^ + 27r/2, where 9 is the rapidity of the particle propagating in the 
loop. 

(b) Vacuum tunneling between two adjacent vacua (/c' = A; ± 1), the loop is the world line 
of a soliton propagating around the cylinder. 

the Hilbert space of the massless free boson introducing an upper energy cutoff (called the 
truncation), which gives a finite dimensional space in finite volume L. The Hamiltonian 
(2.2) can then be represented as a finite numerical matrix and the energy levels evaluated 
by numerical diagonalization. Truncation introduces a source of systematic errors, known 
as truncation errors, which are analogous to discretization errors in lattice field theory. 
Since the wave-functions of the the interacting theory are expanded in terms of lowest- 
lying states of the limiting conformal model (obtained as the high energy limit L — )■ 0), 
truncation errors grow with volume and they are also larger for higher lying levels. In most 
of the calculations in this paper they are too small to be considered, with the exception of 
the issue of resonances treated in section 5. 

The extension of TCSA to the sine-Gordon model was developed in [21] and has found 
numerous applications since then. The Hilbert space can be split by the eigenvalues of 
the topological charge Q (or winding number) and the spatial momentum P, where the 
eigenvalues of the latter are of the form 



In sectors with vanishing topological charge, we can make use of the symmetry of the 
Hamiltonian under 

C: $(a;,t) ^ -$(x,t) 

which is equivalent to conjugation of the topological charge. The truncated space can be 
split into C-even and C-odd subspaces that have roughly equal dimensions, which speeds 
up the diagonalization of the Hamiltonian. 

Using relation (A. 6) we can express all energy levels and matrix elements in units 
of (appropriate powers of) the soliton mass M, and we also introduce the dimensionless 
volume variable / = ML. The general procedure is the same as in [12, 13]: the particle 
content of energy levels can be identified by matching the numerical TCSA spectrum 
against the predictions of the Bethe-Yang equations (3.1). After identification, one can 
compare the appropriate matrix elements to the theoretical values given by (3.5). 

A detailed study of breather form factors based on the formalism explained in subsection 
3.1 has been already performed in [16], therefore it is not repeated here. For the sake of 
concentrating on the essence of the matter, the effect of yU-terms is demonstrated on selected 
examples. In the process of the computation, the analysis was performed for a large number 
of different matrix elements and also for several values of coupling constant ^, so the data 
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shown are just a small, but representative subset of results. In view of the next section 
(dealing with resonances) the most interesting operator to study is 

O = e*i*(°) 

to which therefore all the form factor plots correspond. The numerical data presented 
below correspond to the coupling ^ = 50/311, for which it is instructive to compare the 
scales involved. The masses of the first three breathers in terms of the soliton mass M are 

mi = 0.49973 ■ ■ ■ x M 
ms = 0.96775 ■ ■ ■ x M 
ms = 1.37439 ■ ■ • x M 

The /i-term scales are 





= 0.12486- ■ 


■ X 


M 




= 0.24181 ■■ 


■ X 


M 




= 0.96828 ■■ 


■ X 


M 


2 


= 0.87514 ■■ 


■ X 


M 




= 0.72647 ■■ 


■ X 


M 



Therefore it can be seen that the leading exponential corrections are those driven by the 
splittings B2 —7- BiBi and — t- B1B2, which are exactly the ones described by the 
considerations in subsection 3.2. It is also meaningful to beyond the leading asymptotics 
to evaluate them (i.e. to consider the full solution of the associated Bethe-Yang equations), 
because even 4/i^;^ and 2/x^2 smaller than the next scale mi (although the latter only 
slightly for this particular value of the coupling). 

4.1 Finite volume mass corrections 

One can start by checking whether the mass corrections predicted by eqns. (3.11) for 
B2 and (3.14) for B^ states agree with the numerical energy levels extracted from TCSA. 
Note that there are two copies of zero-momentum B2 and B^ levels in a 2-folded model, 
so the relevant question is whether the finite volume mass gaps evaluated from both of 
them agrees with a single /i-term correction, at least for volumes that are large enough for 
higher exponential terms to be negligible. Indeed, Fig. 4.1 demonstrates that it is so. The 
dashed line shows the infinite volume mass value, while the continuous one is the /i-term 
corrected behavior. The discrete dots show the two mass gaps measured in TCSA. 

Note that the correction is rather small (of the order of at most a few percent) for 
ML > 10. Since the smallest particle mass for the coupling considered is 

mi = 2Msin— ^ 0.5M 
2 

this corresponds to volumes larger than five times the correlation length of the model. 
According to the considerations in [2], this should be enough to neglect these corrections 
in practical computations (in QCD the role of mi is played by the pion mass m^). 

4.2 Vacuum one-particle matrix elements 

For these matrix elements, we plot the dimensionless quantity 

ML) = M-^^/^{0\O\{0})k,L (4.1) 
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Es{L)/M 
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1.1 - 




(a) B2 (b) B, 

Figure 4.2: yU-term corrections to the vacuum-one-particle matrix element at ,^ = 50/311. 
In each case, the dashed line shows the constant behavior expected when neglecting expo- 
nential corrections, the continuous line is the one including the /i-term corrections, while 
the TCSA matrix elements are represented by the two sets of discrete data points. 

where 

IGtt 2(1+0 

is the exact scaling dimension of the operator O. Using formula (3.5) this is naively 
expected to be constant 

ML) = M-^'Ff (0) + O (e-'^^) (4.2) 

where 

(^) = {o\o\B,ie)) 

is the appropriate infinite volume form factor (which is eventually independent of 6 due to 
Lorentz invariance). As demonstrated by Fig. 4.2, however, the exponential corrections 
predicted by (3.15) and the first equation in (3.18), coming from the bound state structure 
of B2 and i?3 lead to observable deviations from this behavior, which agrees quite well with 
the numerical results. Still, one could say that in the regime ML > 10 these are expected 
to be less than 10%, and their presence can be detected from the fact that fk{L) is not 
constant. Note that for each k, there are two one-particle Bk states, and with the two 
vacua this gives four possible matrix elements. However, two of these are very small due 
to suppression by tunneling. 
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Figure 4.3: /i-term corrections to the one-particle-one-particle matrix element at ^ 
50/311 

4.3 One-particle— one-particle matrix elements 

Turning now to matrix elements of the form 



= M-^^/^^/^ i({0}|O|{0})fc,L (4.3) 
we encounter a surprise. According to (3.5), such matrix elements are expected to equal 

/i|fc(L) = M-^FS,itn, 0) + O (e-^^) (4.4) 

where 

FS,{e,,9,) = {0\O\B^i9^)B,i9,)) 
which under the crossing property (3.4) entails that 

/i|fc(L) = M--{B,mO\Bkm + O (e-'^^) 

so this function must be approximately equal to a constant. However, it is very obvious 
from Fig. 4.3 that the exponential corrections are huge (of the order of 100% for k = 2) 
and the constant value (showed by the dashed line) is approached very slowly. It is also 
apparent that the //-term corrected predictions (the continuous lines) from (3.16) and (3.19) 
describe the numerical data very well. For the case of B^ the "naive" //-term prediction 
from the second equation in (3.18) is also shown with dotted line, demonstrating that it is 
very important to get the disconnected parts of the /i-terms right to achieve agreement. 



4.4 Two-particle— one-particle matrix elements 

Let us now study 

/n|.(L) = M~^^Jpu{e, -e)y^ n({-^, ^}|C^|{0})m 
where 6 solves the Bethe-Yang equation 



(4.5) 



miLsinh^ + 5n(2^) = - , 5n(^) = -2 log(-Sii(^)) 
The expectation from the naive finite-volume formula (3.5) is that 



(4.6) 
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/n|2(L) /ii|3W 




Figure 4.4: /i-term corrections to the two-particle-one-particle matrix element at ,^ = 
50/311. The observable deviation for L > 15 in graph (a) is entirely due to numerical 
problems: the raw datum determined from TCSA is the matrix element ({ — ^, |}|C^|{0})2,L; 
which is so small in this range of the volume that it is of the same magnitude as the 
truncation error itself. 



where 

F,'l,{e,,e2,es) = {o\o\B^{e,)B,{e2)Bk{es)) 

Once again, Fig. 4.4 demonstrates that the exponential corrections play an important role, 
and yield corrections of up to 20-30%. This is very important because using the crossing 
property (3.4) one can see that the matrix element measured by the functions fii\k is 

F^,,{t7r-e,t'K + e,o) = {B,{e)B^{-e)\o\BM) 

which for = 3 drives the decay -B3 BiBi, under a perturbation by the operator 

- (C - O^) = sin 

Here we made use of the charge conjugation symmetry 

C: ^{x,t) ^ -^{x,t) 
of the model which ensures that 

{B^{e)B^i-9)\0^\Bsi9s)) = -{B^i9)B^i-e)\0\Bsi9s)) 

because 

COC~^ = for O = e'2* 

and both Bi and B^ are odd under C. According to the data in Fig. 4.4, by omitting the 
/x-term contribution we would significantly underestimate the true infinite volume matrix 
element (the dashed line) on the basis of the TCSA data. Also note that although the 
decay of /i-terms is exponential in the asymptotic regime, in the volume range available in 
the numerical method it could be so slow as to be practically unobservable as displayed in 
Fig. 4.4. What is most surprising is the marked contrast to the case of the energy level 
corrections in Fig. 4.1, which are driven by the same scales fif^ and /i52- 

It is shown in the next section that together with another source of error (also related 
to the /i-terms) this explains the dominant source of systematic error encountered in the 
case of the B^ resonance in the two-frequency sine-Gordon model [11]. 
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5 Decay widths of resonances in the two-frequency sine- 
Gordon model 



Now we turn to an application of the results in the previous section. Consider the two- 
frequency sine-Gordon model defined by the Hamiltonian 



Hdsg = j dx : ]^{dt^f + \{dM : - j dx (^\ : cos/3$ : +C : sin^$ :^ (5.1) 

where normal-ordering is understood treating the field $ as a free massless boson, i.e. 
the theory is treated as a perturbed conformal field theory similarly to (2.2). One can 
consider this model to be a perturbation of sine-Gordon field theory by the integrability- 
breaking coupling C, and using the exact form factors in the C = 0, a so-called form factor 
perturbation theory (FFPT) can be developed [22, 23]. The spectrum and phase structure 
of this model has been studied extensively in the literature [24, 25, 26, 27, 28, 23]. The 
integrability breaking interaction can be assigned the dimensionless parameter 

t = -4^ (5.2) 

This model can be considered as a toy model for hadrons: the C = theory is in our setting 
analogous to QCD, with the breathers being the mesons, which in this limit are all stable. 
Switching on ( makes the mesons unstable, and the ones that are over the threshold 2mi 
are allowed to decay. Consider the decay of B^, which in the rest frame has the kinematics 

with the rapidities of the outgoing BiS given by energy conservation. 

, ^ mg sin^ 

2cosh(^e) = — = — ^ 
mi sm ^ 

Using form factor perturbation theory, the decay width can be calculated [11] 

(rrizY mi I m-j, \ _i 
Km) M \I \ 2mi J ^ 

where the (dimensionless) decay matrix element is given by 

S3ii(^c) = ^-^^ where /3ii(e,) = \Fl^{m,e,,-e,)\ = \Fg^{m,e,,-e, 



where 



^ = sin |$(0) = — (C - O^) with O = e'^ 



In the previous work [11], the prediction (5.3) was compared to numerical TCSA data. 
We developed two methods, of which we only consider the so-called "improved mini- 
Hamiltonian" (imH) method, because it is effectively a direct numerical determination 
of the decay matrix element S3ii(6'c) of the perturbing Hamiltonian. 

The results are shown in Table 5.1 where we use the following parameter for the sine- 
Gordon coupling: 
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0.03022±5 ■ 10-5 


0.0224 


25.946 



Table 5.1: Measured values of S311 using the "improved mini-Hamiltonian" method, com- 
pared to the theoretical prediction labeled by FFPT. The last column contains the volume 
corresponding to the level crossing at zero perturbing coupling. The uncertainties shown 
are crude estimates of truncation errors from TCSA. 



^ V47r ^ 1 



13 ^ 2i?2 _ 1 

(the value ^ = 50/311 used in section 4 corresponds to i? = 1.9). 

Notice that the agreement is not particularly good. The errors are estimates of the 
truncation errors coming from the energy cutoff introduced in TCSA, which is analogous 
to the lattice spacing in QCD. For R > 1.7 this clearly does not explain the deviation 
between the numerical and theoretical results. It was already noted in the original paper 
[11] and it was (correctly) thought to originate from exponential finite size corrections, 
however at that time the tools to test this hypothesis had not yet been developed. 

In subsection 4.4 we already considered this particular matrix element in finite volume 
and noted that the /i-terms did explain the deviation between the naive expectation (3.5) 
and the TCSA data. To investigate whether this helps in the case of the resonance width 
determination, we have to correct errors coming from two sources: 

1. The position of the crossings between the B-^ and the BiBi levels is not where one 
naively expects it to be. Neglecting /i-terms, this would be at the position where 

En(Lo)=m3 (5.4) 

with 

Eii{L) = 2mi cosh^ 

where 6 solves 

miLsinhe + 6n{0) = 2711 (5.5) 
with 1 = 1/2 for the first two-particle level. 

When incorporating the leading /i-term correction to the B^ level, the equation to 
solve is 

EuiL) = Es{L) (5.6) 



where 

with u solving (3.12) i.e. 



E3{L) = mi(l + cos2n) 

gimiLsmh{±m)^^^(^.^)^^^(^2m) = 1 
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Table 5.2: Positions of line-crossings: comparing the measured values and the theoretical 
results obtained without/with the B3 yU-corrections. 
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/i-imH 


S31l{0l) 


1.7 


0.806(4) 


0.868 


1.9 


0.401(1) 


0.401 


2.2 


0.138 


0.134 
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0.0528 


0.0509 
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0.0392 


0.0378 


2.7 


0.0294 


0.0285 



Table 5.3: Comparing the TCSA results for S311 using the "improved mini-Hamiltonian" 
method with /x-term corrections, compared to the theoretical prediction taken at the mod- 
ified kinematical point 61 instead of dc- 

The results are displayed in Table 5.2: note that the yU-term does explain the 
position of the level crossing for i? > 1.9 very well. The remaining discrepancy 
results from further exponential corrections, both to the and the -Bi-Bi level. 
Note that as a result of the difference between (5.4) and (5.6), the rapidity 6 of the 
Bi particles does not equal to 6c, so that means we are not measuring the matrix 
element at the proper kinematical point. As a result, the matrix element relevant to 
the comparison is 

S31l(^l) 

where 61 is the solution of (5.5) at the actual position of the line crossing. 

2. The other effect comes from the yU-term correction for the volume dependence of the 
matrix element, which can be determined from the third formula in (3.18). Notice 
that this requires the knowledge of the five-particle form factor function Fmu. This 
can be used to correct the result obtained from the "improved mini-Hamiltonian" 
method. 

The end result is summarized in Table 5.3. Notice that the agreement between theory and 
numerics is improved by an order of magnitude for R > 1.9, to a few percent deviation 
instead of 10-50 % (the almost exact equality at i? = 1.9 is simply a numerical accident). 
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6 Conclusions and outlook 



6.1 Why are the /i-terms so large? 

The first question to address is: how could the /i-terms be so large in the case of one- 
particle-one-particle and two-particle-one-particle matrix elements? They decay exponen- 
tially by the volume, and therefore are expected to decrease very quickly. As a rule of 
thumb volumes larger than five times the correlation length of the theory are expected to 
be large enough to safely neglect them when calculating to an accuracy of a few percent. 
This is borne out by the finite volume mass corrections examined in subsection 4.1. It is 
important to note that the /i-term predictions for the mass corrections read 

E2{L) = 2mi C0SU2 
E^{L) = mi(l + 2 COSM3) 

with 

«2 = Y + 0(e-'^n^) 
= TT^ + 0(e-'^?2^) 

Note that these enter as arguments of cosine functions, which are very smooth and have 
derivatives of order 1. As a result, their contributions are relatively small for /iL > 1, with 
the asymptotic forms 

£;2(L)=m2-(7i\)V?ie"''"^ + ... 
E3(L) = m3-2(7L)VLe-^'^'' + --- 

What about matrix elements then? The great difference in the yU-term corrected finite- 
volume matrix elements (3.15,3.16), (3.18) and (3.19) is that the parameter u also appears 
in them as kinematical argument of form factors. When the parameters u are expanded 
around their asymptotic value, the correction takes the form of the /i-term exponential 
multiplied by the derivative of the form factor^. Form factors, however, can be rapidly 
varying functions of their arguments. Indeed, they have singularities corresponding to 
bound states and even in the absence of these they always have disconnected contribu- 
tions. These disconnected contributions can appear whenever there are particle rapidities 
displaced by ztt, which according to the crossing relation (3.4) corresponds to a matrix 
elements with one or more particles with identical quantum numbers (including momenta) 
in the two states between which the matrix elements of the local operator are taken. 

The data are indeed consistent with this interpretation. As shown on the example of 
vacuum-one-particle matrix elements, /i-term contributions are similar in magnitude to the 
mass corrections. There are countless other data generated in previous studies (and also 
for testing during this work) for matrix elements between vacuum and multi-particle states, 
and the results are always similar: /i-terms appear, but their magnitude is consistent with 
what is expected from the energy levels. They can only contribute much when their scale 
/i is very small, i.e. for volumes //L < 1 (which can still be a headache in some cases, and 
this provided the motivation for the investigation in [14]). 

However, for matrix elements with multi-particle states on both sides one runs a good 
chance of running into problems with the rapid variation of the form factor appearing in 
the /i-term formulae. Sometimes it is even possible to hit a disconnected piece head-on, 

^We have not written these expansions expHcitly; the interested reader can consult [14] for examples. 
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as happened for the stationary Bi-B^ matrix element for which the theoretical prediction 
(3.19) did include an additional term which introduced an additional factor of L into the 
volume dependence. 

6.2 Implications for the determination of matrix elements 

When matrix elements are to be determined numerically, whether in lattice QCD or in the 
TCSA approach, it is necessary to model the volume dependence of the matrix elements 
in order to extract the infinite volume result one is interested in. The usual way to model 
it accounts for the difference in the normalization between the infinite volume and finite 
volume matrix elements, taken into account by the state density factors in eqn. (3.5). It 
is important to include the corrections by the scattering between the particles in finite 
volume (taken into account by the Bethe-Yang equations (3.1)), because these are of order 
of some inverse power of L. Indeed, this is what Lellouch and Liischer do in their work [2] . 
Exponential finite size corrections, on the other hand, correspond to radiative corrections 
specific to finite volume where the loop winds around because of the periodic boundary 
conditions, and are often neglected. 

In some cases (notably those involving only stationary particles) the properly normal- 
ized matrix elements are predicted to be constant up to exponential finite size corrections. 
So the presence of such corrections can be observed by the fact that the matrix element is 
not a constant, as in subsections 4.2 and 4.3. Therefore one could at least have the idea 
whether there are sizable exponential corrections, and with a bit of luck one may even try 
to fit it and extrapolate to infinite volume. However, in cases involving moving particles as 
in subsection 4.4, this does not work: there is no way of telling the qualitative difference 
between the volume dependence of the dashed and continuous lines in Fig. 4.4 without 
having a theoretical understanding of the /i-terms. 

If one tries to determine the infinite volume form factor without accounting for the 
/i-term, one can easily commit a (systematic) error of the order of 20-50% (for most of 
the parameter space in the examples of two-particle-one-particle form factors considered 
here), and sometimes even errors of the order of 100% (for the Bi — B2 case)! 

Therefore the question is whether it is possible to have a theoretical understanding 
of /i-terms. In the case of integrable models the job is made easier by the existence of 
exact solutions to the form factor bootstrap, in terms of which one can express the /i-term 
contributions by analytically continuing the description (3.6). Note that for the case of 
BiBi — i?3 matrix elements, which is of interest to resonances in the two-frequency sine- 
Gordon model, this requires knowledge of a five-particle form factor, a quite nontrivial 
object. Therefore this does not seem a viable option in the case of theories like QCD, 
where the task is to determine weak interaction matrix elements between hadron states. 

We remark that F-term type exponential finite size corrections are essentially pertur- 
bative when considered in terms of low energy effective field theory; this is already clear 
from the derivation of finite volume mass corrections by Liischer [19]. In lattice QCD, 
they are handled by using chiral perturbation theory to construct the loop integrals, which 
are then put into finite volume using the finite volume perturbation theory formalism 
[29, 30, 31, 32, 33]. The unknown matrix elements entering them are parametrized by 
some low-energy constants, which are to be determined by fitting to volume dependence 
numerically obtained on the lattice. 
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6.3 Implications for resonances 

For resonances, there is not even the silver hning that existed for matrix elements involving 
only stationary particles, where one could at least observe if the volume dependence was 
not quite right. The reason is that their properties can be measured by observing level 
splittings, which occur at particular values of the volume. Even in the two-dimensional 
setting, finding two or more level crossings (more properly these are level avoidances) of the 
same type (meaning the same particle content for the levels) at substantially different values 
of the volume, and with the numerics accurate enough, is a tough problem [11]. However, 
to see whether there are substantial exponential finite size corrections to the decay rates 
one needs to do exactly that: the most straightforward way out is to find level crossings 
at several different values of the volume, and perform an infinite volume extrapolation, or 
at least get an estimate for the systematic error caused by the exponential corrections. 

6.4 Closing remarks 

Without having a handle on the size of the exponential finite size corrections one cannot 
tell whether the observed discrepancy e.g. between a lattice result and an experimentally 
measured decay width is a result of systematic error or a disagreement between theoretical 
model and physical reality. Therefore it is important to have a good theoretical description, 
or at least an estimate of these corrections. One special feature of /i-terms is that in contrast 
to all other exponential corrections which are bounded by the value of the mass gap of the 
theory, /i-terms can decay substantially slower with the volume in the case of weakly bound 
states. Moreover, even when they are expected to decay fast on the basis of the magnitude 
of their characteristic scale /i, their contributions to matrix elements can still be enhanced. 
Measuring their contribution to energy levels is not enough to estimate the errors in matrix 
element and resonance width calculations. 

It is important to stress that the ingredients essential to the observations made in this 
papers are not restricted to the two-dimensional setting. The description of particles as 
composites in subsection 3.2 is just the well-known quantum mechanical treatment of bound 
states as analytic continuation of scattering states to complex values of the momenta. 
In fact, the expression of /i-term corrections to finite volume masses are not that much 
different in higher space-time dimensions [19], the only difference is the appearance of 
some (negative) powers of L (and numerical factors) in front of the exponential correction 
term. Concerning form factors, the matrix elements in quantum field theories of four (or 
indeed any) space-time dimensions share the essential feature with their two-dimensional 
counterparts that they have singularities corresponding to disconnected pieces and bound 
state poles^. Therefore they too can be expected to show rapid variation in kinematical 
domains close to these singularities. In fact, the difficult aspect of modeling the /i-terms in 
higher dimensions is to have the necessary information about many-particle form factors, 
which is so readily provided by the bootstrap in integrable models. 
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A Sine- Gordon breather form factors 



A.l Spectrum and S matrix 

The fundamental excitations of the sine-Gordon model (2.1) are a doublet of soliton/antisoliton 
of mass M. Their exact S matrix can be written as [5] 



where 



and 
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Soi9,0- 



sinh ( I 



sinh 
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2 sm I I 



sinh 
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(A.l) 



2 sinh 2|* cosh f 

(the latter representation is valid for any value of n G N and makes the integral represen- 
tation converge faster and further away from the real 9 axis). 

Besides the solitons, the spectrum of theory contains also breathers B^, with masses 



rrir 



2M sin ■ 



(A.2) 



The breather-soliton and breather-breather 5*- matrices are also explicitly known [5]. 



A.2 Breather form factors 

We only consider exponentials of the bosonic field $. To obtain matrix elements containing 
the first breather, one can analytically continue the elementary form factors of sinh-Gordon 
theory obtained in [35] to imaginary values of the couplings. The result is 

FiV.i(^i,...,^n) = (0|e^'^^*(°)|5i(ei)...i?i(^„)> 
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and 
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(A.4) 



gives the minimal BiBi form factor^, while cr^"^ denotes the elementary symmetric poly- 
nomial of n variables and order k defined by 
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(A.5) 



is the exact vacuum expectation value of the exponential field [36], with M denoting the 
soliton mass related to the coupling A in (2.2) via [37] 
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Formula (A. 3) also coincides with the result given in [38]. Form factors containing higher 
breathers can be obtained using that is a bound state of Bi and -B„_i; therefore 
sequentially fusing n adjacent first breathers gives Following the lines of reasoning 
of Appendix A of the paper [11] one obtains 
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(A.8) 



■^The formula for the function v is in fact independent of N; choosing N large extends the width of the 
strip where the integral converges and also speeds up convergence. 
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is the BiBk — )> -Bfc+i coupling, defined as the residue of the appropriate scattering amph- 
tude: 



i7r(fc + l)g 
2 



sinh 9 + i sin ^ ^ ' sinh + z sin ^ 
sinh 6* — z sin ^^"^ ^ sinh 6' — z sin ^ 
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